Radio Identification of Sub-mm Sources in the Hubble Deep Field 
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ABSTRACT 

Determination of the epoch dependent star-formation rate of field galaxies is 
one of the principal goals of modern observational cosmology. Recently, Hughes et 
al. (1998) using the SCUBA instrument on the James Clerk Maxwell Telescope, 
report the detection of a new population of heavily dust enshrouded, star-forming 
galaxies at high redshifts (z >2), dramatically altering the picture of galaxy evo- 
lution. However, we show that this interpretation must be treated with caution 
because of ambiguities in the identification of the host galaxies. Based on our 
deep, high resolution 1.4 GHz obervations of the Hubble Deep Field, we suggest 
alternate identifications to the sub-mm detections. These identifications argue 
for a substantially lower redshift to the sub-mm population with a consequential 
lowering of the z >2 sub-mm/far infrared luminosity density and global star- 
formation rate. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
starburst, radio continuum 
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1. Introduction 

The data used to make a quantitative estimate of the global star-formation history typi- 
cally consist of optical multi-color imaging combined with comprehensive spectroscopy (Lilly 
et al. 1996). These data are then used to infer the distances of individual galaxies and (with 
the aid of stellar synthesis models) determine their intrinsic spectral energy distributions. In 
this manner active star-formation rates can be calculated (Madau et al. 1996, Connelly et 
al. 1997). The primary limitation of this technique is its sensitivity to dust obscuration, the 
effects of which in evolving galaxy populations are both uncertain and controversial (Madau 
et al. 1996, Calzetti 1997, Meuer et al. 1997). 

An alternative method for the study of star-forming galaxies is to observe directly the 
reprocessed UV light emitted by the dust. Depending on the dust temperature (30 K - 60 
K), this emission peaks between 50-100 fxm (restframe) and in many galaxies comprises the 
bulk of the bolometric luminosity. Additionally, the thermal far infrared radiation (FIR) 
is not subject to further obscuration, so uncertain extinction corrections are avoided. At 
high redshifts, the FIR emission will be shifted into the sub-mm band. A consequence of the 
steep Rayleigh- Jeans tail of the thermal dust emission coupled with a positive dust emissivity 
index, a (where the emissivity, e oc v a ) is to bring higher flux into the sub-mm band so the 
cosmo logical dimming effect is offset. Thus a starburst galaxy of a given luminosity should 
have essentially the same observed flux density between z = 1 and 10 (Blain et al. 1993). 

Closely related to the FIR emission in starburst galaxies is the radio continuum. In 
normal galaxies (i.e., without a powerful AGN), the centimeter radio luminosity is dominated 
by diffuse synchrotron emission believed to be produced by relativistic electrons accelerated 
in supernovae remnants. At shorter wavelengths free-free emission from HII regions may 
contribute substantially. Although the radio emission is linked to active star-formation by 
different physical mechanisms than that of the FIR, there is a tight correlation between the 
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FIR and radio luminosity of a starburst (Helou et al. 1985, Condon et al. 1991). Radio 
observations are only sensitive to recent starburst activity in a galaxy (and the formation 
of its O and B stellar populations) since the thermal and synchrotron radiation dissipate 
on physical time scales of 10 7 — 10 s yr. In this sense the radio luminosity of a starburst is 
a true measure of the instantaneous rate of star-formation in a galaxy, uncontaminated by 
older stellar populations. Because galaxies and the inter-galactic medium are transparant 
at centimeter wavelengths, radio emission is a sensitive measure of star-formation in distant 
galaxies. The current deep Very Large Array (VLA) radio surveys with sensitivites of a few 
microjansky are capable of detecting star-forming galaxies to z ~1.5 (Richards et al. 1998). 

An 850 ^m survey of the HDF with the SCUBA detector on the James Clerk Maxwell 
Telescope (JCMT) detected five sources in a confusion limited image (Hughes et al. 1998) 
with 15" resolution. Tentative optical identifications are all with putative starbursts at z > 1 
and with star-formation rates (SFR) of 400-1200 M & yr~ x (we assume h = 0.5, q Q = 0.5) 
In this letter we compare our deep radio images of the Hubble Deep Field with the SCUBA 
images and suggest alternate optical counterparts to the sub-mm sources. 

2. Radio Observations 

A 5.4' (FWHM) field containing the HDF has been imaged at 8.5 GHz with the VLA 
with an rms sensitivity of 1.8 //Jy. The observing technique and data reduction are discussed 
in Richards et al. (1998). We collected 40 additional hours of data at 8.5 GHz in June 1997. 
The new combined image has a rms sensitivity near the field center of about 1.6 ^Jy with a 
resolution of 3.5". 

During November 1996, we obtained 42 hours of VLA data in its A-array on the HDF 
at 1.4 GHz. The subsequent 1.4 GHz image of the HDF has an effective resolution of 
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1.9" and an rms sensitivity of 7.5 //Jy with a firm detection limit of 40 /iJy. A total 
of 381 radio sources at 1.4 GHz have been catalogued within 20' of the HDF and are re- 
ported on elsewhere (Richards 1998). Within the Hubble Deep Field, there are nine radio 
sources detected in complete samples at either 1.4 GHz and/or 8.5 GHz, while an additional 
seven 8.5 GHz sources comprise a supplementary sample as described by Richards et al. 
(1998). The 8.5 and 1.4 GHz images of the HDF are available at the world-wide web site: 
www.cv.nrao.edu/ jkempner/vla-hdf/. 



3. Association of Sub-mm and Radio Sources 

We inspected our radio images to determine if any of the SCUBA sources have possible 
radio counterparts. Our first step was to align the radio and SCUBA position frames. The 
VLA coordinate grid is within 0.1" of the J2000/FK5 reference frame at both 8.5 and 1.4 
GHz (Richards et al. 1998). In order to tie the JCMT coordinate grid to this system, we 
have assumed that the radio and sub-mm sources are associated with the same galaxy, as 
discussed in the introduction. 

The relative rms positional errors for the sub-mm sources should be of order 1-2" (based 
on the 15" SCUBA beam size and the signal to noise ratio of individual detections), howver 
the uncertain effects of source confusion in the sub-mm images likely makes this an under- 
estimate. In addition, the absolute registration of the SCUBA image is unknown a priori, 
although Hughes et al. (1998) quote a value of 0.5", while Smail et al. (1998) report typ- 
ical values of 3" for their SCUBA images. Thus we chose to search for any radio object 
either in the 1.4 GHz complete sample or in the 8.5 GHz catalog of Richards et al. (1998) 
within a 10" error circle around each of the sub-mm source positions. Possible radio associ- 
ations were apparent for HDF850.1 (3651+1226 and 3651+1221), HDF850.2 (3656+1207), 
and HDF850.4 (3649+1313). Two of these radio sources were detected at both 1.4 and 8.5 
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GHz as part of independent and complete samples (3651+1221 and 3649+1313). These two 
radio sources in particular also have high signifigance ISO 15 fim counterparts (Aussel et al. 
1998) indicating these systems may contain substantial amounts of dust and hence be lumi- 
nous FIR galaxies. Based on the association of 3651+1221 and 3649+1313 with HDF850.1 
and HDF850.4, respectively, we suggest a SCUBA coordinate frame shift of 4.8" west and 
3.8" south. 

The position shifts, within 1.2" in the (VLA+ISO) vs. SCUBA positions for both 
HDF850.1 and HDF850.4, suggests that the VLA/ISO/SCUBA source alignment is not 
accidental. However this large registration offset of ~6" is much greater than the 0.5" reg- 
istraion accuracy quoted in Hughes et al. (1998). Either the radio/ISO emission is not 
associated withe the same galaxies as the SCUBA sources, or the SCUBA observations have 
a large registration error. 

4. Optical Identification of Radio/Sub-mm Sources 

Since the radio source positions are much more accurate (0.1-0.2") than the sub-mm 
positions and as the HDF contains a high surface density of optical objects (typically 20 
per SCUBA beam), we now use the radio data to make the secure identifications with 
optical counterparts. Table 1 presents plausible radio counterparts to the sub-mm sources of 
Hughes et al. (1998). The first line gives the position of the SCUBA source after translation 
to the radio coordinate frame with plausible radio counterparts and their suggested optical 
identification given in following lines. 

HDF850.1 : We present in Figure 1, the 1.4 GHz overlay of the Hubble Deep Field 
centered on the shifted SCUBA position. The precise 1.4 GHz radio position suggests the 
optical identification is with the faint low-surface brightness object to the immediate north 
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of the brighter foreground disk system. Richards et al. (1998) suggested that this optical 
feature might be associated with the z = 0.299 galaxy 3-659.1. However, we note the presence 
of a separate low surface brightness galaxy at z = 1.72 (3-633.1; Fernandez-Soto et al. 1998) 
located approximately 2" to the northwest. 

If radio source 3651+1221 is the most obscured part of a larger galaxy 3-633.1 at z 
= 1.72 (Fernandez-Soto et al. 1998), the implied radio luminosity (L L4 = 10 25 W/Hz) is 
substantially higher than that of the most extreme local starbursts (e.g., ARP 220) and 
suggests that this object may contain an AGN. 

HDF850.2 : Based on the optical/radio positional coincidence, this 3.5 a radio source 
has a 90% reliability of being associated with a 1=23.7 distorted galaxy (Barger et al. 1998), 
according to the analysis of Richards et al. (1998). We identify the 850 /im detection with 
this galaxy of unknown redshift. The UGR band photometry (Hogg 1997) on this galaxy 
suggests that it is likely at z < 3. Figure 2 shows the 1.4 GHz overlay of the optical field. 

HDF850.3 : The radio data does little to clarify the identification of the sub-mm source 
in this field. Figure 3 shows a 4 a radio source located 4" from the position of the 850 /im 
detection. However, before the shift the SCUBA source position is in good agreement with 
the position of the bright disk system 1-34.0, which is also included in the supplemental ISO 
catalog of Aussel et al. (1998). This galaxy also has a weak 3a radio detection. The 0.485 
redshift of this object implies a star-formation rate of 80 M© yr" 1 from the radio luminosity 
(Salpeter initial mass function integrated over 0.1-100 M Q ), although the presence of an 
AGN cannot be ruled out. At present the data cannot discriminate between these two 
possible radio/sub- mm associations. 

HDF850.4 : The radio source 3649+1313 is associated with the spiral galaxy 2-264.1 at 
a redshift of 0.475. ISO sources from the complete catalogs of Aussel et al. (1998) and Rowan- 
Robinson et al. (1997) have been assoicated with this radio source. This galaxy is likely 
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part of the larger structure at 0.475 which contains 16 galaxies with spectroscopic redshifts 
(Cohen et al. 1996). At least one of these galaxies (2-264.2) lies at a small projected distance 
(less than 30 kpc) and suggests dynamic interactions may be triggering the radio/sub-mm 
activity. Although the SCUBA detection may be a blend of emission from several galaxies in 
this crowded field (see Figure 4), the radio emission is clearly confined to the central galaxy 
2-264.1. Richards et al. (1998) estimate a SFR = 150 M Q yr" 1 . 

Can the SCUBA source HDF850.4 plausibly be associated instead with the HDF optical 
galaxy 2-399.0 as claimed by Hughes et al.l If we take those authors' estimate of the FIR 
luminosity logi L 6 o Mm , = 12.47 L & for this galaxy, the FIR-radio relation (Condon et al. 
1991) predicts an observed 1.4 GHz flux density of about 300 //Jy, clearly in excess of our 
upper limit of 23 fiJy (3 a). We find the identification of SCUBA source HDF850.4 with 
HDF galaxy 2-399 to be dubious and instead identify HDF850.4 with 2-264.1. 

HDF850.5 : There is no 1.4 GHz radio emission apparent to the 2 a limit of 15 /iJy in 
this field. Optically there are only two plausible identifications for the 850 /an source, HDF 
galaxies 2-395.0 and 2-349.0 (see Figure 5). The redshift of 2-395.0 is 0.72 (Fernandez- Soto 
et al. 1998). The radio flux limit on this galaxy can exclude a SFR > 60 M & /yr. The 
other possible identification (2-349.0) is at an unknown redshift. The non-detection of this 
sub-mm source at radio or optical wavelengths coupled with the fact that this is the weakest 
source in the Hughes et al. (1998) catalog suggests this source may be spurious. We also 
note that this sub-mm source is located only 12" from sub-mm detection HDF850.4 and 
hence may suffer from confusion. 
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5. Conclusions 

Of the five 850 /jm sources in the HDF, two are solidly detected at radio wavelengths, 
while two are probable detections. Two of these identifications are possibly with z ~ 0.5 
starbursts (HDF850.3 and HDF850.4). The other two detections (HDF850.1 and HDF850.2) 
must have redshifts less than 1.5 or be contaminated by AGN based on radio luminosity 
arguments. This radio analysis suggests the claim, based on sub-mm observations alone, 
that the optical surveys underestimate the z > 2 global star-formation rate are premature. 
On the other hand, the z < 1 star-formation history may have been underestimated if a 
significant fraction of the sub-mm population lies at relatively low redshift. 

In the absence of high resolution sub-mm imaging capability, it is necessary to rely on 
plausible radio counterparts of sub-mm sources in order to provide the astrometric accuracy 
needed to make the proper optical identifications. Only complete redshift samples of the 
sub-mm population coupled with diagnostic spectroscopy and high resolution radio data will 
allow for calculation of the epoch dependent sub-mm luminosity function and its implication 
for the star-formation history of the universe. 

We thank our collaborators K. Kellermann, E. Fomalont, B. Partridge, R. Windhorst, 
and D. Haarsma for a critical reading of an earlier version of this work. We appreciate useful 
conversations with J. Condon and A. Barger. Support for part of this work was provided 
by NASA through grant AR-6337.*-96A from the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS5-2655, and by NSF grant AST 93-20049. 
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Figure Captions 



Fig. I. — The greyscale shows a 20" x 20" HDF I-band image (Williams et al. 1996) contain- 
ing the SCUBA detection HDF850.1. The contours correspond to 1.4 GHz emission at the 
-2, 2, 4 and 6 a level (a = 7.5 /iJy). The three sigma position error circle for HDF850.1 is 
shown after shifting to the VLA coordinate frame. The original position of HDF850.1 taken 
from Hughes et al. (1998) is denoted by the diamond. The ISO detection is marked with 
a cross with three sigma position errors (Aussel et al. 1998). The radio emission is clearly 
confined to the optical feature north of the bright spiral. Radio source 3651 + 1221 may be 
the most obscured part of a larger galaxy 3-633.1 at z — 1.72 (Fernandez-Soto et al. 1998). 

Fig. 2. — The greyscale shows a 20" x 20" ground-based I-band image taken from Barger 
et al. (1998) centered on the position of SCUBA source HDF850.2. The 1.4 GHz radio 
contours are drawn at -2, 2, 3 and 5 a. The symbols are the same as for Figure 1. We 
identify HDF850.2 with the 3.5 o radio source 3657+1159. 

Fig. 3. — The greyscale corresponds to optical I-band emission in the field of SCUBA source 
HDF850.3 (20" on a side) as taken from the HDF. Radio contours at 1.4 GHz are drawn at 
the -2, 2 and 4 a level. Intriguingly, there is a 4.2 a radio 'source' located 4" from HDF850.3. 
The probability of this being a chance coincidence is 20% based on the surface density of 
4(7 radio sources in the field. If HDF850.3 is associated with the radio source then this is a 
blank field object to Iab = 27 (object lies in the less sensitive PC). However, an ISO source 
from the supplemental catalog of Aussel et al. (1998) is also in the field and associated with 
the bright disk galaxy 1-34.0 (Williams et al. 1996). It is difficult to discrimante between 
these two possible sub-mm identifications with the present data. The symbols are the same 
as for Figure 1. 

Fig. 4. — Radio 1.4 GHz contours drawn at the -2, 2, 4 and 6 a level are overlaid on the 
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HDF I-band image centered on the position of HDF850.4 (20" on a side). A 15 yum detection 
from the complete catalog of Aussel et al. (1998) has been associated with this radio source 
and suggests likely starburst activity in the disk galaxy. The symbols are the same as for 
Figure 1. 

Fig. 5. — Radio 1.4 GHz contours drawn at the -2 and 2 a level are overlaid on the HDF 
I-band image in the field of HDF850.5 (20" on a side). This is the one sub-mm source in the 
HDF which has no plausible radio counterpart. The symbols are the same as Figure 1. 



This figure "richards_figl.gif" is available in "gif" format from: 
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Table 1: Radio Identifications of Sub-mm Sources 



Name RA 1 DEC 1 HDF ID 2 S 8 5G Hz 3 S 14G Hz 3 Ub 4 z 5 

(AtJy) (^Jy) (mag) 

(1) (2) (3) (4) (5) (6) (7) 



HDF850.1 
3651+1221 



12:36:51.63±0.10 
12:36:51.653±0.021 



62:12:22. 5±0. 7 
62:12:21.41±0.15 



3-633.1? 



17 



49 



25.9 



;i.72)? 



HDF850.2 

3655+1159 



12:36:55.99+0.20 
12:36:55.859+0.033 



62:12:00.0+1.4 
62:12:01.57+0.23 



<6 



26 



23.9 



HDF850.3 

3644+1258 
3644+1304 



12:36:44.06+0.21 
12:36:44.627+0.028 
12:36:44.612+0.036 



62:12:59.5+1.5 
62:12:58.15+0.19 
62:13:04.73+0.25 



1-34.0 



<5 
<5 



32 
23 



>27 
21.3 



0.485 



HDF850.4 
3649+1313 



12:36:49.68+0.23 
12:36:49.656+0.019 



62:13:12.1+1.6 
62:13:13.16+0.13 



2-264.1 



14 



49 



20.6 0.475 



HDF850.5 



12:36:51.29+0.25 



62:13:15.4+5.4 



1 The SCUBA positions are given after translation to the VLA coordinate frame. 

2 The most probable optical identification of the radio source in the HDF catalog (Williams 
et al. 1996) based on positional coincidece. 

3 The integrated flux density corrected for the attenuation of the VLA primary beam. In a 
few places upperlimits are given. 

4 The integrated I magnitude of the identification as taken from the HDF catalog (Williams 
et al. 1996). We have transformed to the AB system using the convention (HDF magnitude - AB) 
= 0.8. 

5 The spectroscopic redshift (Cohen et al. 1996, Lowenthal et al. 1997) of the identification. 
Values within parentheses are photometric estimates of the redshift (Fernandez-Soto et al. 1998). 



